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Abstract 

A concise method has been formulated for identifying a set of forces needed to constrain 
the behavior of a mechanical system, modeled as a set of particles and rigid bodies, when 
it is subject to motion constraints described by nonholonomic equations that are inherently 
nonlinear in velocity. An expression in vector form is obtained for each force; a direction is 
determined, together with the point of application. This result is a consequence of expressing 
constraint equations in terms of dot products of vectors rather than in the usual way, which 
is entirely in terms of scalars and matrices. The constraint forces in vector form are used 
together with two new analytical approaches for deriving equations governing motion of a 
system subject to such constraints. If constraint forces are of interest they can be brought into 
evidence in explicit dynamical equations by employing the well-known nonholonomic partial 
velocities associated with Kane’s method; if they are not of interest, equations can be formed 
instead with the aid of vectors introduced here as nonholonomic partial accelerations. When 
the analyst requires only the latter, smaller set of equations, they can be formed directly; 
it is not necessary to expend the labor to form the former, larger set first and subsequently 
perform matrix multiplications. 
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1. Introduction 


Motion constraints imposed on a mechanical system are described with nonholonomic 
(nonintegrable) constraint equations, whereas configuration constraints are expressed with 
holonomic constraint equations. Two examples of motion constraints with which the reader 
may be familiar are the condition of rolling, which is the absence of slipping, and the restric- 
tion on velocity imposed by a sharp-edged blade. These constraints are sometimes described 
with equations written in the matrix form au + f3 = 0, where u is a column matrix of motion 
variables u\, . . . ,u n . Motion variables, also referred to as generalized speeds, are in general 
linear combinations of the time derivatives of generalized coordinates, q\, ... ,q n . The distin- 
guishing feature of such equations is that they are linear in the motion variables. However, 
one may consider motion constraints that must be described by relationships that are inher- 
ently nonlinear in the motion variables, having the form f(qi, . . . ,q n ,Ui, . . . ,u n ,t ) = 0. In 
Ref. [1] Bajodah et al. review some of the literature dealing with nonlinear nonholonomic 
constraint equations and consider it important to study them because they can arise in 
connection with servo-constraints or program constraints when a control system enters the 
picture. As explained in Refs. [2] and [3], such constraints are enforced by application of 
control forces as opposed to the forces present when bodies and particles come into contact 
with one another, as is the case with classical, passive constraints. 

Golubev states in Ref. [4] that, as of yet, there is no example of a passive mechanical 
device that can compel a motion constraint described by an equation that is nonlinear in 
velocity. Roberson and Schwertassek note in Ref. [5] that all known motion constraints im- 
posed on purely mechanical systems can be expressed with relationships that are linear in 
velocity variables. Unfortunately, the relationships in such situations are often artificially 
teased into nonlinear forms to create contrived examples used to illustrate a proposed pro- 
cedure. For instance, a nonlinear equation is devised in Ref. [6] to describe the constraint 
imposed on a rolling disk. The well-known Appell-Hamel mechanism is studied and dis- 
cussed, for example, in Refs. [1] and [7] - [12]. It is recognized in Refs. [1] and [8] - [12] 
that the constraints imposed on this mechanical system can be expressed with linear rela- 
tionships, but despite this the mechanism is used in Refs. [1], [11], and [12] to demonstrate 
the application of methods for dealing with nonlinear nonholonomic constraint equations. 
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In Refs. [13] and [14], Zekovich offers several examples of passive mechanical systems in 
which the constraints are described with nonlinear nonholonomic constraint equations. In 
what follows it is shown that the associated constraints can in fact be expressed with linear 
nonholonomic equations. Another example studied in Refs. [8], [15], [16], and [17] involves 
a device proposed by Benenti in Ref. [18]. However, a purely mechanical system is involved 
and therefore, according to the observations in Refs. [4] and [5], the nonlinear nonholonomic 
equation used to describe the constraint must be regarded as contrived. 

Whenever a motion constraint can be expressed entirely with linear nonholonomic con- 
straint equations, it should be dealt with accordingly. Any number of approaches can be 
used to deal with the equations in their linear form; the exercise of cajoling such equations 
into a nonlinear appearance serves no useful purpose. The new approaches contained in this 
paper, and the examples of their application, are concerned strictly with inherently nonlinear 
nonholonomic constraint equations. 

The literature contains several instances of motion constraints described by nonholonomic 
equations that are inherently nonlinear in velocity. Perhaps the simplest case, provided by 
Golubev in Ref. [4], involves a single particle P that is subject to a uniform gravitational 
field and moves in a vertical plane fixed in an inertial reference frame N. The magnitude of 
the velocity N v p of P in N is to remain constant. The particle thus constrained serves as a 
model of a robot manipulator tip used to spray-paint a wall or polish a surface. Variations of 
this problem are studied in Refs. [19]— [21] . A familiar example proposed by Appcll, in which 
P moves in three dimensions, is discussed in Refs. [7], [15], [19], [22], and [23]. Special cases 
of Appcll’s problem are examined in Refs. [20] and [24], Control of an inverted pendulum 
constitutes an example studied in Refs. [15] and [16]. A thin rigid rod moves in a vertical 
plane in the presence of a uniform gravitational field, with the lower end of the rod always 
in contact with a horizontal line. The system is referred to as Marie’s servomechanism; as 
proposed in Ref. [7], an actuator controls the horizontal displacement of the rod’s lower end 
according to some control law in order to keep the rod vertical. An earlier paper by Huston 
and Passerello (Ref. [25]) considers the more general case of balancing a pole whose lower 
end remains in contact with a horizontal plane, while the pole is otherwise free to move in 
the space above the horizontal plane. 

The forthcoming developments in this paper are carried out for the most part in terms 
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of vectors. These quantities are used also in expressing the main results, and discussing the 
contributions of the work. By vector we mean a basis-independent quantity having direction 
and magnitude, such as position, velocity, acceleration, or force, involved in the application 
of elementary principles of dynamics to study motion taking place in three-dimensional 
space. Other examples of a vector include partial velocities and partial angular velocities 
associated with advanced principles of dynamics. We do not mean a row or column matrix 
whose elements consist of three basis-dependent scalar measure numbers of a vector. Nor do 
we have in mind a matrix containing more than three scalar elements, such as a collection of 
generalized forces, or a row or column matrix considered from the viewpoint of linear algebra 
to belong to an rt- dimensional tangent space, orthogonal space, etc. 

In Ref. [26], a comprehensive, consistent, and concise method is established for identifying 
a set of forces needed to constrain the behavior of a mechanical system modeled as a set of 
particles and rigid bodies. The method is exercised in Ref. [27] with an example involving a 
configuration constraint, and a motion constraint expressed with an equation that is linear 
in velocity. The purpose of this paper is to apply the method to constraints described by 
nonholonomic equations that are inherently nonlinear in velocity. (It is to be understood 
that the term “velocity,” used in the general case of a system of particles, subsumes “angular 
velocity” in the special case in which a subset of particles makes up a rigid body. The 
term “acceleration” likewise encompasses an angular counterpart.) An essential feature of 
the method consists of expressing constraint equations in vector form rather than entirely 
in terms of scalars and matrices as is customary. A constraint equation that has been 
differentiated once or twice with respect to time, so that it contains the acceleration of a 
point or the angular acceleration of a rigid body, is said to be written at the acceleration 
level. Likewise, a constraint equation at the velocity level is one that has been differentiated 
at most once, so that it contains the velocity of a point or the angular velocity of a rigid 
body. It so happens that the method discussed in Refs. [26] and [27] can be applied whenever 
constraints can be described at the acceleration level by a set of independent equations that 
are linear in acceleration; therefore, it is applicable to constraint equations that are nonlinear 
in velocity when written at the velocity level. 

The method in question yields expressions in vector form for constraint forces, and for 
torques of couples formed by constraint forces (hereafter referred to as constraint forces and 
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torques). Thus, the directions of these vectors are identified, together with the specific point 
at which a constraint force must be applied, and the particular body upon which a constraint 
torque must be exerted. Such information about the vector quantities is of interest in its 
own right, and is to be preferred over the information contained in a matrix whose elements 
are scalar generalized constraint forces. In the process of constructing generalized constraint 
forces, information about the direction, magnitude, and point or body of application of 
constraint forces and torques becomes lost; in principle, each generalized constraint force is 
a sum of contributions from every constraint force and torque acting on a mechanical system. 
Although generalized constraint forces can be computed in a straightforward manner from 
knowledge of constraint forces and torques, usually it is impractical to invert the process 
and recover the original information about constraint forces and torques from generalized 
constraint forces. 

Anderson is concerned in Ref. [28] with configuration constraints and with motion con- 
straints described by nonholonomic equations that are linear in the motion variables. Al- 
though such constraints are not the direct subject of the present investigation, Anderson 
makes an observation that is nevertheless relevant to our discussion. Often, a Lagrange 
multiplier or undetermined multiplier used to treat a constrained system is not related in a 
clear way to any particular constraint force or torque. In the method introduced here, each 
multiplier has a straightforward relationship to a constraint force or torque. 

The emphasis in this paper is on analytic derivation of equations of motion that do or do 
not contain evidence of forces and torques needed to impose motion constraints described 
with inherently nonlinear nonholonomic equations. This stands in contrast to methods of 
computational dynamics, where the object is numerical formulation and solution of equa- 
tions of motion. With knowledge of constraint forces and torques obtained by inspection of 
constraint equations written in vector form, and the two new approaches developed here, the 
analyst can form explicit equations of motion by hand or with the aid of symbolic algebra 
software. Equations that do not contain evidence of constraint forces can be formed directly; 
they need not be obtained from numerical manipulations of equations in which evidence of 
constraint forces is present. 

The remainder of the paper is organized as follows. First, a treatment of nonlinear 
nonholonomic constraint equations is undertaken in Sec. 2 for a generic system of particles; 
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the results are applicable whether or not a subset of particles makes up a rigid body. The 
method of Ref. [26] is used to identify directions of constraint forces and the particles to which 
they must be applied. The constraint forces are used together with extensions to Kane’s 
method (Ref. [30]) to obtain two new ways of deriving dynamical equations of motion. The 
first of these is useful when one is interested in the time histories of the constraint forces; 
it produces dynamical equations that contain evidence of the constraint forces needed to 
satisfy the nonlinear nonholonomic constraint equations. On the other hand, the second 
approach can be used when one is not interested in the constraint forces but requires explicit 
dynamical equations governing the motion of the constrained system; constraint forces are 
not in evidence in the minimal equations of motion obtained with this approach. The 
novelty in the second case rests in the use of nonholonomic partial accelerations rather than 
the nonholonomic partial velocities employed in Kane’s method. The methods proposed 
in Sec. 2 are first compared in Sec. 3 to two of the approaches in the existing literature, 
and then applied in Sec. 4 to an example in which the velocities of two particles must 
remain perpendicular. The resulting equations of motion are solved numerically. Constraint 
forces are identified in Sec. 5 for two other examples in which the velocities of two particles 
must either remain parallel, or equal in magnitude. In connection with Appell’s particle, 
a constraint force is identified in Sec. 6; a second demonstration of the two approaches for 
obtaining equations of motion is performed, and the equations are compared to existing 
results. Finally, Sec. 7 contains the essential steps that must be taken to extend the ideas 
presented in Sec. 2 from a discussion in terms of a system of particles to the practical case 
in which a subset of the particles makes up a rigid body. Concluding remarks are supplied 
in Sec. 8. 

2. Equations of Motion for Complex Nonholonomic Systems 

Thomas R. Kane has been developing and extending an approach to solving problems 
in dynamics for the past five decades. An early paper, Ref. [29], and two highly influential 
books, Refs. [30] and [31], are but three of the notable publications by Kane and his col- 
leagues. Kane’s method enjoys widespread application in the areas of multibody dynamics 
(Ref. [32]), dynamics of complex spacecraft (Refs. [31] and [33]), and robotic devices (Ref. 
[34]). Two well-known computer programs, discussed in Refs. [35], [36], and [37], employ 
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computer algebra in carrying out Kane’s method to derive equations of motion specific to 
the system of interest, and subsequently to create software for numerical solution of the 
equations. 

Kane’s method can be used to construct dynamical equations of motion belonging to a set 
that is minimal in number; in other words, there are as many equations as there are degrees 
of freedom in the mechanical system. In pursuing this approach, constraint forces may be 
treated in the same way as other contact forces and distance forces, or the constraint forces 
may be left out of the picture, because in either case they do not contribute to the equations 
of motion. For this reason such forces are referred to as noncontributing. On the other hand 
it is important to note that, when time histories of these forces are of interest, Kane’s method 
contains provisions for bringing them into evidence selectively. In this case the dynamical 
equations are greater in number than the degrees of freedom of the system. Whether or 
not the constraint forces contribute depends on whether dependent and independent motion 
variables, or only independent motion variables, are included in expressions for velocities of 
points, and angular velocities of rigid bodies, when such expressions are inspected to identify 
vectors known as partial velocities and partial angular velocities. 

Kane’s method is set forth in full detail in Ref. [30]. The analyst chooses to form one of 
the following three sets of dynamical equations of motion, depending upon whether a system 
S is subject to configuration constraints and/or motion constraints, and what constraint 
forces, if any, are of interest. 


f + + 

= 0 

(r = 1,.. 

. , n + M) 

(1) 

F r +f; 

= 0 

(r = 1,.. 

■ , n) 

(2) 

F r +K 

= 0 

(r = 1,.. 

. ,n — m) 

(3) 


The quantities without a superscript * are referred to as generalized active forces, whereas 
the quantities with a superscript * are known as generalized inertia forces. Configuration 
constraints imposed on S are described by M independent holonomic constraint equations, 
and motion constraints are represented by m independent nonholonomic constraint equations 
that are linear in the motion variables. The configuration of S in a Newtonian reference 
frame N is described by n generalized coordinates. 

The superscript + is not used in Ref. [30] ; we employ it here to signify equations formed 
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according to Secs. 4.9, 6.3, and 7.6 therein for the purpose of bringing into evidence forces 
associated with holonomic constraint equations. Such forces make contributions to F r + but 
not to F r or F r . Constraint forces that must be applied to satisfy any linear nonholonomic 
equations contribute to F r and to F r + , but not to F r . The apparatus of Ref. [30] deals only 
with holonomic constraint equations and linear nonholonomic constraint equations; there 
are no provisions whatsoever for dealing with nonlinear nonholonomic constraint equations. 
Consequently, F r as dealt with in Ref. [30] contain no evidence of constraint forces of any 
kind. In what follows we propose extending Kane’s method in order to deal with nonlinear 
nonholonomic constraint equations. 

It is instructive to recall that configuration constraints are, in general, expressed at 
the position level with nonlinear holonomic constraint equations. However, when these 
relationships are expressed at the velocity level they are linear in the velocity vectors or, what 
is the same, linear in the motion variables as shown in Ref. [26]. Similarly, motion constraints 
in general are described at the velocity level by nonlinear nonholonomic constraint equations 
but, when expressed at the acceleration level, they are linear in the acceleration vectors. 
In other words, when written in scalar form the latter relationships are linear in the time 
derivatives of motion variables. 

It is also important to remember that the partial velocities used to form Eqs. (1) are 
obtained from velocity expressions that do not account for configuration constraints, whereas 
the partial velocities employed in constructing Eqs. (2) and (3) are collected from velocity 
expressions that do account for configuration constraints. It is precisely for this reason that 
forces associated with holonomic constraint equations contribute to F r + but not to F r or F r . 

Two important conclusions follow from these observations. First, because inherently 
nonlinear nonholonomic constraint equations written at the acceleration level are linear in 
acceleration vectors, the forces needed to satisfy those constraints can be identified with 
the approach described in Ref. [26]. Second, those forces can be brought into, or left out 
of, evidence in equations of motion by making use of partial accelerations obtained from 
acceleration expressions that respectively do not, or do, account for the associated motion 
constraints. 

Suppose that a simple nonholonomic system S (Ref. [30]) is made up of particles P ±, . . . , P v . 
The configuration of S in a Newtonian reference frame N is described by generalized coor- 


dinates qi,...,q n , and the motion of S is characterized by independent motion variables 
ui, ... ,u p , where p = n — m. Suppose further that S is subject to i independent nonlinear 
nonholonomic constraint equations 

h,( N v p ‘,..., N v p ",t) = 0 (s = l,...,«) (4) 

where N v Pi is the velocity of particle Pj ( i = 1, . . . , u) in N, and where t denotes time. 
In this case S is referred to as a complex nonholonomic system. Differentiation of these 
relationships with respect to t in N yields 

^2 N & Pi • W is + Z s = 0 (s = 1,...,£) (5) 

i= 1 

where W is are vector functions of q 1 , . . . , q n , ui , . . . , u p and t in N, and Z s are scalar functions 
of the same variables. The acceleration of Pi in N is represented by N a Pi . When these 
independent relationships are satisfied the motion variable time derivatives ui, ... ,u p are no 
longer independent of one another, as discussed presently. 

By virtue of Newton’s second law, Eqs. (5) have certain implications regarding the con- 
straint forces needed for their satisfaction. According to Ref. [26] one can inspect Eqs. (5) 
and conclude that constraint forces are given by 

C is = X s W is (i = l,...,z/; s = 1,...,£) (6) 

where X s are scalar multipliers whose time histories may, or may not, be of interest. The 
constraint force C i S is evidently parallel to W !S , and in general it must be applied to Pi in 
order to satisfy the constraint equations (5). One may use terminology from Golubev’s Ref. 
[4] to refer to C j S as an ideal servoconstraint force. A non-ideal servoconstraint force could be 
formed as AW is + Cj_, where C_l is a force perpendicular to W is . When available actuators 
are incapable of exerting an ideal servoconstraint force, it may be possible to satisfy the 
constraint with a non-ideal force. In this paper we limit the discussion to the ideal case. 
The fundamental definition of Kane’s generalized active forces involves the dot product of 
two vectors; C i S is one such vector. Knowledge of the direction and point of application of 
C i S , and its relationship to A s , is important for its own sake. It is at least as important as 
having a collection of generalized constraint forces in hand, if not more so. The technique 
of inspecting Eqs. (5) systematically establishes the direction and point of application of a 
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constraint force very soon after a constraint equation is available at the acceleration level in 
vector form, generally much sooner and with less labor than when working with constraint 
equations written entirely in terms of scalars and matrices. 

The constraint forces formed according to Eqs. (5) and (6) do make contributions to F r , 
therefore Eqs. (3) can be used if such constraint forces are of interest. A new set of equations 
are developed for use in place of Eqs. (3) in the event that these constraint forces are not of 
interest. 

The contributions of Cj S to F r are revealed by expressing Eqs. (3) in terms of fundamental 
definitions from Ref. [30], 

Fr + K = J2 Y’ • (Ri - mi N a p ‘) 

2=1 

= j^ N v r Pi -('f i + J^A s W is -m i N a Pi )=0 (r = l,...,p) (7) 

i= 1 V .9=1 / 

where F r , F*, and N v Pi respectively denote the rth nonholonomic generalized active force 
for S in N, nonholonomic generalized inertia force for S in N, and nonholonomic partial 
velocity of P ?; in N. The mass of P, is indicated by m, . The resultant R.j of all contact forces 
and distance forces acting on P t is regarded as the sum of the constraint forces X) s=1 A S W i S 
that must be applied to ensure satisfaction of Eqs. (5), added to the resultant of all other 
forces, fj. Equations (7) together with Eqs. (5) furnish the number of relationships needed 
to solve for the unknown quantities hi, . . . , u p , Ai, . . . , A^. One employs these relationships if 
the time histories of Ai, . . . , Ag are of interest. 

A reduced or minimal set of dynamical equations to which Cj S do not contribute is given 
by 

F r + F; = jr N a r p '- (Vj + ^ A s W is — rrii N a Pi 

i= 1 V s=l 

= jr N a Pi • (fj - mi N a Pi ) =0 (r = 1, . . . , c) (8) 

Z — 1 

where 

c = p — £ (9) 

is the number of degrees of freedom of S in N. When speaking of F r and F* it is convenient 
to refer to them, respectively, as the rth nonholonomic generalized active force for S in N 
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and the rth nonholonomic generalized inertia force for S in N, but the double tilde notation 
should be used to indicate they have been formed with N a Pi , the rth nonholonomic partial 
acceleration of P* in N, rather than ;V v r f * . When one is not interested in time histories 
of Ai, . . . , A £, one can construct Eqs. (8) directly by forming the dot products indicated in 
the second line rather than the first line. Moreover, Eqs. (7) need not be constructed first. 
Directly forming Eqs. (8) thus eliminates the need for assembling a constraint Jacobian 
and an orthogonal complement, and subsequently using the latter matrix to annihilate the 
former. 

It is important to realize that the nonholonomic partial accelerations N a Pi in Eqs. (8) 
are distinct from the nonholonomic partial velocities N v Pi in Eqs. (7). In addition, one must 
have practical instructions for obtaining the vectors N a Pi . Finally, it is essential to point 
out that use of the vectors N a p ‘ makes it unnecessary, in general, to include the constraint 
forces in Eqs. (8). 

The acceleration of P t in N can be written uniquely in terms of «!,... , u p , 

p 


N _ 


= E 


N a Pi u r 


+ 


N Pi 


(i = 1, 


v) 


( 10 ) 


r= 1 


and, also uniquely, in terms of the independent motion variable time derivatives 


N a Pi - 


= E 


N~Pi 


U r 4" 


N~Pi 


(f = 1, . . . , Z/) 


( 11 ) 


r = 1 


Equations (10) and (11) are analogous to Eqs. (2.14.2) and (2.14.4) in Ref. [30], where it 
is established that a holonomic partial velocity v r is distinct from a nonholonomic partial 


velocity v . Similarly, the partial acceleration 




is decidedly different from the nonholo- 


nomic partial acceleration 


Nx Pi 


because the right hand member of Eqs. (11) involves only 


the independent motion variable time derivatives. 

Equations (10) can be obtained from Eq. (2.14.4) of Ref. [30] by differentiation with 
respect to t in N, in which case the partial acceleration N a Pi is seen to be identical to the 
nonholonomic partial velocity of Pj in N, 


JV„ Pi A N~ Pi 


a r % = 


(i — 1, . . . ,i/; r — 1, ... ,p) 


and the acceleration remainder N a t p ‘ is defined to be 




P /N j 

A P / d 


N 


H 


= £ 7 


Nx r Pi 


d 


NT r P l 


r — 1 


dt 


Ur + j t v * 


(• = i, ■■■,<') 


( 12 ) 


(13) 
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Substitution from Eqs. (10) into (5) gives 

E (E 'V 1 • w is ) u r + E "V ■ w ?;s + = 0 (s = 1 , . . . ,£) (14) 

r = 1 \t=l ) i= 1 

The coefficients of h r and the remaining terms can be abbreviated respectively by means of 
two definitions, 

a sr = E (s = 1, ... ,^; r = 1, ... ,p) (15) 

i = 1 

and 

% = + E "V • W ? ; s (s = l,...,£) (16) 

2=1 

where a sr and are functions of qi, . . . ,q n ,Ui, . . . ,u p , and the time t. These definitions 
allow Eqs. (14) to be rewritten in a form that is linear in the time derivatives of the motion 
variables 

p 

EawMr + 7s = 0 (s = l,...,£) (17) 

r = 1 

These relationships express the dependence of l time derivatives of the motion variables, 
say u c+ 1 , . . . , iip, on the remaining ones hi, ... ,ii c . ft is assumed that these independent 
equations can in fact be solved for h c+1 , . . . , u p in terms of hi, . . . , ii c . The dependent motion 
variable time derivatives are written in terms of the independent ones in a manner analogous 
to Eqs. (2.13.1) of Ref. [30], 

C 

hc+r ^ ] A r sU s T B r (r 1, . . . , £) (1^) 

S=1 

With a relationship for N a Pi in hand having the form of Eqs. (10), one simply embeds 
the acceleration level constraint equations by rewriting h c+ i, . . . ,u p in terms of hi, . . . ,h c 
to obtain an expression in the form of Eqs. (11). Nonholonomic partial accelerations N dL r p ‘ 
are subsequently obtained in the same way as partial velocities, namely by inspecting the 
resulting relationship for acceleration to determine the vector coefficients of h r for r = 
1, . . . ,c. 

When dealing with simple nonholonomic systems and the associated constraint equations 
(2.13.1) of Ref. [30], the analyst chooses which p of u ±, . . . , u n to regard as independent; of 
course, the remaining motion variables are then regarded as dependent. The choice is made 
during the process of deriving explicit equations of motion. The same is true in the case 
of Eqs. (18) here; the analyst chooses which c of hi, . . . ,h p are considered independent. In 
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neither case is the decision based on the result of numerical procedures used in connection 
with the computational method of coordinate partitioning discussed in Refs. [38] and [39]. 
One does not, for instance, “take advantage of the numerical structure of the Jacobian 
matrix” (Ref. [39]). As Anderson notes in Ref. [28], coordinate partitioning is an iterative, 
computationally expensive procedure that cannot be used in explicit symbolic formulation 
of equations of motion. 

The remainder of this section is devoted to a discussion of the contributions of the 
constraint forces C is (i — 1, . . . , zq s — 1, . . . , €) to generalized active forces. The contribu- 
tions to F r (r = 1, . . . ,p) are examined first, and consideration of the contributions to F r 
(r — 1, ... ,c) follow. 

Nonholonomic generalized active forces for S in N, F r , are defined by Eqs. (4.4.1) in Ref. 
[30] as the sum of dot products of pairs of vectors: 

Fr = jZ N v r Pi • R* (r = 1, . . . ,p) (19) 

2=1 

Let C i represent the resultant of the constraint forces C,;. s applied to Pi in order to ensure 
satisfaction of Eqs. (5), so that 

e £ 

C; E c - (i = (20) 

S=1 S=1 

The resultant of all contact forces and distance forces acting on Pi can then be regarded 
as the sum of the constraint force, Cj, and the resultant of all other forces, f ? . Hence, F r 
is made up of contributions (F r )c from the constraint forces acting on S and (F r )jr from all 
other forces acting on S, 

Fr = (• Fr) C + (F r )T = E ^ • Q + £ ^ {r = l, ... , P ) (21) 

2=1 2=1 

The contribution from the constraint forces can be singled out, and it is given by 

{F r )c = £ N v r Pi • E = E (r — 1, ... ,p) (22) 

2=1 S = 1 5=1 

where, keeping in mind Eqs. (12), a sr has the same meaning as in Eqs. (15). As is true when 
obtaining any generalized active force by using the techniques of Ref. [30] , the recommended 
approach is to form the dot products indicated in Eqs. (19), (21), and (22). 
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It can be shown that, in general, the constraint forces C make no contribution to any 
of F r . A general proof is omitted in the interest of brevity; however, the result can be stated 


as 



t 


i 


= o (r = 1, . . . , c) 


(23) 


Therefore, it can be concluded that 



(24) 


i = 1 


In words, Eqs. (23) state that motion constraints described by inherently nonlinear nonholo- 
nomic constraint equations require the application of forces that make no contributions to 
any of the nonholonomic generalized active forces F r . The utility of this result in practice is 
that, when directly forming Eqs. (8) for a particular system, the constraint forces C, may 
be included in R, or they may be omitted; in either case they will not contribute to F r . 

3. Comparison with Other Methods 

The works cited in the reference list present many ways of constructing equations of mo- 
tion for systems subject to motion constraints described by nonlinear nonholonomic equa- 
tions. We compare the methods proposed in Sec. 2 with two of those existing approaches. 

In Refs. [40] and [41], Udwadia and Kalaba describe a method for obtaining general 
equations of motion for discrete mechanical systems. The dynamical equations are referred 
to variously as nonminimal, unreduced, or full order, because there are more equations than 
there are system degrees of freedom. Another distinguishing feature of the equations is 
that they are free of Lagrange multipliers or any other unknowns representing the constraint 
forces. A critical step in the derivation of the nonminimal equations is observed to be the use 
of constraint equations that have been differentiated an appropriate number of times so that 
they are expressed at the acceleration level; the result is that holonomic and nonholonomic 
systems are treated in a unified way. The nonholonomic constraint equations expressed 
at the velocity level can be either linear or nonlinear in the time derivatives of generalized 
coordinates. The constraint equations expressed at the acceleration level need not be linearly 
independent. 
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The equations of Udwadia and Kalaba are expressed in terms of q r and q r rather than the 
more general quantities u r and u r . Therefore, in the interest of making a comparison, we in- 
voke the limitation that u r is chosen to be simply q r . As noted earlier, Eqs. (7) together with 
Eqs. (5) [or their alternative forms, Eqs. (17) or (18)] can be solved for the unknown motion 
variable time derivatives, and the multipliers, in terms of the motion variables, generalized 
coordinates, and time. If one subsequently eliminates the multipliers from Eqs. (7), the 
resulting relationships are equivalent to Udwadia and Kalaba’s nonminimal, multiplier-free 
equations for the case of independent acceleration level constraint equations. A disadvantage 
of Udwadia and Kalaba’s approach, then, is that time histories of the multipliers are un- 
available, even in the event they are of particular interest. To be sure, generalized constraint 
forces are available with their method; however, this state of affairs is unsatisfactory for the 
reasons mentioned in Sec. 1. Even with generalized constraint forces in hand, one is not able 
to separate out the individual vectorial constraint forces and torques acting on the particles 
and rigid bodies of the system, or determine to which particles and bodies the individual 
forces and torques must be applied. Our method for obtaining such useful information, cen- 
tered around Eqs. (5) and (6), contains no counterpart in the work of Udwadia and Kalaba. 
Finally, they offer nothing in the way of a minimal equation set such as our Eqs. (8). That is 
to say, when constraint forces (even generalized constraint forces) are of no interest whatso- 
ever, their approach involves unnecessary overhead in deriving and solving a greater number 
of equations of motion than is absolutely necessary. The number of unnecessary equations 
of motion will in that case be equal to the number of constraint equations. 

Huston and Passerello (Ref. [25]) were the first to approach the matter of extending 
Kane’s method to deal with nonlinear nonholonomic constraint equations; their work is 
refined in Ref. [11]. A similar viewpoint for dealing with linear nonholonomic constraint 
equations is presented in Refs. [42] and [10]. 

There are certain concepts that the exposition in Sec. 2 has in common with that of 
Ref. [11], The authors of that work recognize constraint equations that are nonlinear at the 
velocity level become linear at the acceleration level, and they note the relationship between 
partial acceleration and nonholonomic partial velocity expressed in Eqs. (12). They make 
use of these observations to form equations of motion that are equivalent to Eqs. (7), and 
form generalized constraint forces that are expressed with the final term in Eqs. (22). It 
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is pointed out that the undetermined multipliers can be eliminated and a reduced set of 
equations of motion can be obtained. 

There exist a number of differences between what is presented here and in Ref. [11], 
In that work, the development is restricted to motion variables that are each defined as 
the time derivative of a single generalized coordinate. Remainder terms such as or 

N v t Pi needed to account for prescribed motion are not included in the formulation. The 
development requires partial velocities to be expressed in a vector basis fixed in an inertial 
reference frame, which is not necessarily convenient or efficient. In contrast, the motion 
variables used here are fully general linear combinations as in Eqs. (2.12.1) of Ref. [30], 
velocity remainder terms A v ( P ' are included [see Eqs. (13)], and all partial velocities (for 
that matter, all vectors) introduced herein are considered basis-independent quantities just 
as they are in Ref. [30]. 

In Ref. [11], equations containing the multipliers are formed first; the multipliers are sub- 
sequently eliminated and a reduced set of equations of motion similar to Eqs. (8) is obtained 
by premultiplication with an orthogonal complement matrix. (An analogous approach is 
taken in Refs. [10] and [42] in connection with linear nonholonomic constraint equations.) 
As is well known, an orthogonal complement is not unique. In simple problems an orthogonal 
complement can be obtained analytically, as in Ref. [11], Usually, however, it is produced 
numerically via the zero-eigenvalue theorem, singular value decomposition, QR decompo- 
sition, successive multiplication of Householder transformations, etc. As noted earlier, the 
Appcll-Hamel mechanism is used to illustrate the method proposed in Ref. [11] even though 
it involves contrived nonlinearity in nonholonomic constraint equations. 

The present work puts forth two signihcant advances over the material in Ref. [11], First, 
information about the direction and point of application of constraint forces is obtained by 
inspecting constraint equations written in vector form at the acceleration level. As demon- 
strated in Sec. 7, the direction and body of application of a constraint torque can be obtained 
in the same way. In Ref. [11] the undetermined multipliers are related in a clear way to scalar 
generalized constraint forces, but not to constraint forces and torques in vector form. Sec- 
ond, it is discovered here that nonholonomic partial accelerations can be used to construct 
Eqs. (8) directly and analytically. This approach circumvents the need to form Eqs. (7) first 
and, afterwards, carry out what are usually two numerical procedures, namely production 
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and application of an orthogonal complement. The absence of orthogonal complements is a 
desirable feature common to the methods of Ref. [30] and this work. There is no introduction 
of the nonholonomic partial acceleration in Ref. [11], or of the nonholonomic partial angular 
acceleration that is defined in Sec. 7. In contrast to nonunique orthogonal complements, 
the nonholonomic partial accelerations and nonholonomic partial angular accelerations pro- 
posed here are unique once a set of independent motion variable time derivatives has been 
chosen, and they are formed by the same definite process of inspection used to obtain partial 
velocities and nonholonomic partial velocities. 

4. Two Particles with Perpendicular Velocities 

An example is provided to illustrate application of Eqs. (7) and (8) to form equations of 
motion in which constraint forces respectively are and are not in evidence. A system of two 
individual particles is subject to a requirement that the velocity in a Newtonian reference 
frame N of one particle must remain perpendicular to the velocity in N of the other particle. 
The associated nonholonomic constraint equation is inherently nonlinear. Implementation 
of the constraint would require the sort of computations that are associated with a control 
system, as well as ideal actuators and sensors; thus, the example features a servo-constraint. 
The demonstration is followed by discussion of a similar example from the literature in which 
the constraint is imposed by purely mechanical means, and it is shown that the nonholonomic 
constraint equation can in that case be expressed as a linear relationship. 

Two pucks moving on an air-bearing table fixed in a Newtonian reference frame N are 
modeled as particles P\ with a mass of mi, and P 2 with a mass of m2- Let two orthogonal 
unit vectors fi! and n 2 be fixed in N and define the plane of the table, and let unit vector 
n 3 = hx x n 2 be normal to the plane. An external force fi = oqrii + <r 2 h 2 is applied to Pi 
whereas a force f 2 = 03111 + 04n 2 is applied to P 2 . The motion of this system is regarded 
as unconstrained. Suppose that the velocities iV v P| and ;V v p ' 2 of Pi and P 2 in N are to be 
constrained such that they must remain perpendicular at all times. 

Let mi = 1 kg, m 2 = 2 kg, and let fi and f 2 be characterized by the constants oq = 1.0 
N, cr 2 = 0 N, 03 = 1.0 N, and 04 = 0 N. At t = 0 the velocities of Pi and P 2 in N are given 
by N v Pl = 0.311! + 0.4n 2 m/s, and N v F - = 0.4hi — 0.3n 2 m/s. The initial position vectors 
p, from a point O fixed in N to P t are given by p x = lrq — 2n 2 m, and p 2 = lhi + 2n 2 m. 
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The remainder of this section is divided into four parts. First, a constraint equation 
is written in the form of Eqs. (5) and subsequently inspected to identify constraint forces 
according to Eqs. (6). This exercise can yield important information about the constraint 
forces, even without forming or solving equations of motion according to (7) or (8). Second, 
Eqs. (7) are employed to produce dynamical equations of motion in which the constraint 
forces play a part, and these equations are solved numerically together with kinematical 
differential equations. Third, an alternative set of equations of motion in which constraint 
forces do not play a part, are formed by carrying out the steps indicated in Eqs. (8), and 
results of a numerical solution are discussed. Finally, a very closely related published example 
involving a classical mechanical system is examined. 

4-1. Identification of Constraint Forces 

The constraint can be expressed by the relationship 

V 2 - V‘=0 (25) 

This constraint equation is nonlinear in the velocity vectors because more than one velocity 
appears in a dot product; it is also nonlinear in motion variables, as will become apparent. 
Differentiation with respect to t in N brings the constraint equation to the acceleration level, 
where it is seen to be linear in the acceleration vectors because only one such vector appears 
in each dot product. 

N a P2 • N v Pl + N a Pl • N v P2 = 0 (26) 

With Eqs. (5) and (6) in mind, it can be concluded that the constraint requires application 
of the forces 

C 2 = A jv v P i , C 1 = \ n v P2 (27) 

to P 2 and P\ respectively. A constraint force can be applied to a puck, for example, by four 
orthogonally mounted thrusters. The constraint forces Ci and C 2 need not be of equal mag- 
nitudes because the constraint does not require ;V v p ' 2 and ;V v P| to be equal in magnitude. 
The constraint force C\ is perpendicular to C 2 when the constraint is satisfied. All of this 
valuable information concerning the vectors Ci and C 2 , including their relationship to A, is 
obtained by inspecting Eqs. (26) rather than by attempting to infer it from examination of 
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generalized constraint forces. The vector forms in Eqs. (27) are required for forming the dot 
products indicated in Eqs. (7). 


4-2. Constraint Forces In Evidence 

Equations of motion containing evidence of Ci and C 2 can be derived according to Eqs. 
(7). The unconstrained system possesses four degrees of freedom in N, thus the motion can 
be characterized by four motion variables defined operationally as 

N v Pl = rqn 4 + u 2 h 2 , N v p ' 2 = u 3 rq + n 4 n 2 (28) 


These relationships are inspected to identify the vector coefficients of tq, u 2l u 3 , and w 4 ; that 
is, the nonholonomic partial velocities 


V 1 = nr, 


N v Pl = n 2 , 


n~ Pi _ 


= 0, 


N~Pi 


= 0 


(29) 


v- p 2 _ 


q =0, A v 2 P2 = 0, "vi*=n 1 , N v p ' 2 = h 2 (30) 

The partial velocities are referred to as nonholonomic, and the notation N v Pi is used to 
indicate that the expressions in Eqs. (28) would have accounted for any nonholonomic con- 
straint equations linear in the motion variables, had any such equations been applicable. 
Dynamical equations of motion formed according to Eqs. (7) are readily written as 


miUi — (7\ + A u 3 , miu 2 = o 2 + A m 4 , m 2 ii 3 = cr 3 + Arq, rri 2 u 4 = a 4 + A u 2 (31) 


The constraint equation expressed at the velocity level in vector form by Eq. (25) becomes, 
in scalar form, 

UiU 3 + u 2 Ui = 0 (32) 

This relationship is nonlinear in the motion variables. As pointed out earlier, Ref. [30] 
contains no provisions for dealing with such a constraint equation, therefore it cannot be 
used when forming familiar holonomic partial velocities or nonholonomic partial velocities. 
It is for this reason that nonholonomic partial accelerations are introduced in this paper; 
these vectors can be used to construct equations of motion devoid of A, as is demonstrated 
shortly. Now, the constraint equation at the acceleration level is linear in the time derivatives 
of the motion variables, 

u 3 ii\ + ?i 4 h 2 + U\U 3 + u 2 iii = 0 (33) 
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An analytical solution of the linear system of equations (31) and (33) for the five unknowns 
is manageable, and is given by 


iii 


^ _ mi(o- 3 ni + (T4M2) + m 2 (criM 3 + (J2M4) 
mi{ui 2 + u 2 2 ) + m 2 (u 3 2 + U4 2 ) 

&i + Am 3 . <7 2 + \u± . cr 3 + \u\ , C4 + Am 2 

, U-2 = , U 3 = , U4 = 

mi mi mo m2 


(34) 

(35) 


The configuration of P\ and P 2 in N is described by four generalized coordinates introduced 
operationally as 


Pi = gihi + g 2 n 2 , p 2 = <?3 ni + g 4 n 2 (36) 

Four kinematical differential equations are given simply by 


q r = u r (r = 1, 2, 3, 4) 


(37) 


The dynamical and kinematical differential equations are integrated numerically with a 
variable step-size algorithm, using an absolute error of 1 x 10~ 8 and a relative error of 1 x 10 “ 7 . 
The unconstrained trajectories (A = 0) of Pi and P 2 are displayed in the upper left of Fig. 1, 
to be compared to the constrained trajectories shown in the upper right. It is clear that ;V v P| 
and N v p ' 2 are becoming parallel in the absence of constraint forces, whereas they remain 
perpendicular when Ci and C 2 are applied. A time history of A is shown in the lower left 
of Fig. 1. The constraint requires A v P2 to remain perpendicular to ;V v P| ; hence, the cosine 
of the angle between the two vectors calculated as cos 9 = A v P2 • ,v v Pl /( v v p -|| A v Pl |), 
which should be 0, can be used as a measure of the failure of the numerical solution to satisfy 
the constraint. As seen in the lower right of Fig. 1, the solution meets the constraint very 
well. 


4-3. Constraint Forces Not In Evidence 

Although use of Eqs. (7) has been demonstrated first, one can of course bypass these 
relationships completely in favor of Eqs. (8) if a time history of A is not of interest. One can 
virtually eliminate the small error evident in the time history of cos 9, and obtain dynamical 
equations of motion in which A does not appear, by appealing directly to Eqs. (8). First, 
the accelerations in N of P\ and P 2 are expressed as 

N a Pl = Miiii + u 2 n 2 , N a Pa = fi 3 ni + u 4 h 2 (38) 
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Figure 1: Two Particles with Perpendicular Velocities 


The motion variable time derivatives iti, ii 2 , and u 3 can be chosen as independent. This 
leaves W 4 as dependent, and one then substitutes from Eq. (33) to arrive at 

N a Pl = Mini + u 2 n 2 , N a P ' 2 = u 3 hi — ~(u 3 ui + w 4 w 2 + uiu 3 )h 2 (39) 

u 2 

The nonholonomic partial accelerations of P\ and P 2 in N are identified as the vector coef- 
ficients of hi, ii 2 , and u 3 , 


N a[ 1 = hr, 


% Pl = n 2 , 


NxPi 


N~P 2 _ U:i ^ 

a i — n 2; 

U 2 


U 4 „ 

a 0 ' z = n 2 , 

U 2 


N~P 2 _ 

x 2 


= 0 


N~ Po - M 1 - 

a , 2 = ni n 2 

u 2 


(40) 

(41) 


These vectors are evidently fewer in number than, and distinct from the nonholonomic partial 
velocities in Eqs. (29) and (30). Once they are in hand, nonholonomic generalized active 
forces for S in N can be formed according to the expressions 


Fr = 


Pi 


(f 1 + A Ar v P2 )+ N a r P2 ■ (f 2 + \ N v Pl ) (r = 1,2,3) 


(42) 
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The first of these is given by 


Fi — hx • (fi + \ N v P2 ) - — n 2 • (f 2 + \ N v Pl 

u 2 

— cri + A u 3 - — (cr 4 + Au 2 ) 

U 2 

U 3 

— 01 04 

U 2 


( 43 ) 


Similarly, 

F 2 = cr 2 - — cr 4 (44) 

U 2 

F 3 = CT 3 - — cr 4 (45) 

U 2 

The multiplier A is clearly absent from F] , F 2 , and F 3 , and thus the constraint forces Ci and 
C 2 do not contribute to the equations of motion. Considering the result stated in Sec. 2, 
one would be justified in omitting Ci and C 2 from Eqs. (42) and (43), and thereby reducing 
the labor involved in forming dot products. The nonholonomic generalized active forces are 
obtained without first constructing relationships according to Eqs. (7), and without forming 
a Jacobian matrix or its orthogonal complement, and multiplying them together. 
Nonholonomic generalized inertia forces are given by 


F* = N a Pl 


N 

-mi a 


Pl ) + N a P2 


-m 2 N a Pz 


) (r = 1,2,3) 


(46) 


or 


u 3 


F* = —mini - m 2 — ~(u 3 iii + m 4 m 2 + u\u 3 ) 


u 2 Z 

, ( U 3 \ 2 

mi + m 2 * — 

\u 2 J 


U 3 Ui . 

U\ - m 2 — 71 -U 2 - m 2 - 


UlU 3 . 


u 2 ‘ 


U 2 ‘ 


U 3 


(47) 


«4 


F* = -miii -2 - m 2 — ^(m3«i + « 4 M 2 + “ 1 %) 


u 2- 


U 3 U 4 


= ~m 2 — y u i ~ 
u 2 z 


mi + m 2 


u 2 - m 2 - 


U 1 U 4 


u 2 - 


u 3 


(48) 


Ml 


F* = -m 2 u 3 - m 2 — o (m 3 m 4 + m 4 m 2 + M1M3) 


m 2 


UiU 3 . 


= -m 2 — —ui — m 2 9 
u 2 z u 2 2 


m 4 m 4 

m 2 - m 2 


1 + 1 ^) 2 

u 2 J 


u 3 


(49) 
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The mass matrix associated with these equations of motion is symmetric. After expressing 
u 4 as —U 1 U 3 /U 2 as required by Eq. (32), the dynamical equations of motion F r + F* = 0 
(r = 1,2,3) and the kinematical differential equations (37) are integrated numerically using 
the initial conditions given in the problem statement. The paths of P\ and P2 are identical 
to those shown in the upper right plot of Fig. 1, and the absolute value of cos 9 remains less 
than 7.64 x ICE 17 throughout the simulation. 

4-4- A Classical Mechanical System 

I 11 Refs. [13] and [14] Zekovich provides examples in which velocities of two particles are 
to remain perpendicular to one another. However, an additional configuration constraint is 
imposed on Pi and P 2 ; they are connected by a “fork” that allows relative translation along 
the line joining Pi and P 2 . bi other words, Pi is regarded as fixed in a rigid body B, and 
a prismatic joint makes it possible for P 2 to move on B. A relationship having the form 
of Eq. (32) is given, and put forth as an example of a nonlinear nonholonomic constraint 
equation. However, the nonlinearity is contrived. The development in Ref. [13] is greatly 
simplified by working with a set of motion variables to be defined presently; furthermore, 
they are used to show that the relevant nonholonomic constraint equations can be written 
as linear expressions. 

Let perpendicular unit vectors bi and b 2 be fixed in B such that they lie in the plane of 
motion of Pi and P 2 , and bi is in the direction of the prismatic joint that permits P 2 to slide 
on B. Unit vector b 3 is perpendicular to bi and b 2 , and to the plane of the motion. Four 
motion variables are introduced operationally by writing N v Pl = Wibi + u 2 b 2 , N uj b = U 3 b 3 , 
and B v Pl = u 4 bi. The angular velocity of B in N is denoted by n uj b , and the velocity 
of P 2 in B is indicated by p v Pz . Hence, N v P2 = {ui + u^hi + (w 2 + q 4 u 3 )h 2 , where q 4 
is the distance between Py and P 2 . The perpendicular velocity constraint is expressed as 
jV v P2 • N v Pl = u x (ui + u 4 ) + u 2 {u 2 + q 4 u 3 ) = 0. 

Zekovich begins the analysis by attaching a sharp-edged circular disk, or blade, at Pi with 
the edge perpendicular to bi; the resulting constraint is expressed linearly as A v Pl • bi = 
Ui = 0, and the corresponding Eq. (8) in Ref. [13] is likewise linear. With u\ = 0, the 
velocity constraint is rewritten as A v P2 • = u 2 (u 2 + q 4 u 3 ) = 0, which corresponds to 

Eq. (9) of Ref. [13]. Zekovich then notes the constraint can be satisfied in either of two 
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ways. The first possibility is imposition of the constraint expressed by the linear equation 
jY v Pl • b 2 = U 2 = 0, in which case P\ is fixed in N and the blade at P\ is no longer 
necessary. The second possibility also involves a constraint described by a linear relationship 


perpendicular blades is in keeping with the result of Eqs. (27), although it contradicts the 
direction of R 2 indicated in Fig. 3a of Ref. [13]. 

5. Other Examples 

Other restrictions on the motion of two separate particles give rise to nonholonomic 
constraint equations that are inherently nonlinear. Constraint forces required to ensure that 
the velocities in N of the two particles remain parallel, or equal in magnitude, are discussed 
briefly. This is followed with a mention of two examples involving a single particle. 

First consider the requirement that ;Y v Pl and N v p ' 2 be parallel to each other. Allow 
the plane containing ,Y v Pl and A v p ‘ 2 to be oriented arbitrarily in iV ; without loss of three- 
dimensional generality, define the unit vector h 3 to be perpendicular to this plane. The 
constraint can then be expressed as follows. The vector n 3 x ,Y v P| is perpendicular to n 3 
and to ,Y v P] by construction; therefore, requiring ,Y v p ' 2 to be parallel to ,Y v P| is the same 
as requiring 


This constraint equation is observed to be nonlinear in the velocity vectors because more 
than one velocity appears in a dot product. Differentiation with respect to t in N brings the 
constraint equation to the acceleration level, where it is seen to be linear in the acceleration 
vectors. 


A v P2 • b 2 = u -2 + q^u 3 = 0; such a restriction can be imposed by fixing a blade at P 2 
with the edge orthogonal to b 2 . The presence of perpendicular constraint forces exerted by 



(50) 



(51) 


In view of Eqs. (5) and (6), the constraint requires application of the forces 


C 2 = A(n 3 x N v Pl ), Ci = — A(n 3 x N v P2 ) 


(52) 


to P 2 and Pi respectively. The constraint forces Ci and C 2 need not be of equal magnitudes 
because the constraint does not require N v p ' 2 and A v Pl to be equal in magnitude. Moreover, 
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Ci and C 2 may have the same direction or opposite directions depending on whether the 
directions of ;V v P| and A 'v Pl are opposite or the same. As is the case in the example in 
Sec. 4, important information about constraint forces is obtained by inspecting a constraint 
equation written at the acceleration level in vector form. Extracting the same information 
from generalized constraint forces would be significantly more arduous. The relationship 
between the multiplier and the two constraint forces is clear-cut. 

The first example in Refs. [13] and [14] is similar to the preceding situation, but an 
additional configuration constraint is imposed on P 1 and P 2 ; they are connected by a rod 
of fixed length 2 L. It is said that the requirement of parallel velocities can be achieved 
in practice by attaching at the rod’s midpoint a blade that is perpendicular to the rod. 
A relationship is given with the form of Eq. (50) written entirely in terms of scalars, and 
offered as an example of a nonlinear nonholonomic constraint equation. However, in this 
instance the nonlinearity is contrived because the constraint dictated by the blade can in 
fact be described by a linear nonholonomic constraint equation. There appears to be some 
recognition of this in Ref. [13]. The directions of the constraint forces obtained in Eqs. (52) 
are seen to be the same as those indicated in the diagram on the right side of Fig. 2a in Ref. 

[13]. 

Next, suppose that N v Pl and N v p ' 2 are required to have equal magnitudes rather than 
parallel directions or perpendicular directions. The constraint can be expressed by the 
relationship 

N v P2 • N v P2 - N v Pl ■ N v Pl = 0 (53) 

which is nonlinear in the velocity vectors. At the acceleration level, the constraint equation 
is linear in the acceleration vectors, 

N a P2 . N v P 2 _N a P i . N vPl = q ( 54 ) 

According to Eqs. (5) and (6), the constraint requires application of the forces 

C 2 = X N v P2 , C 1 = -\ N v Pl (55) 

to P 2 and Pi respectively. It is seen that Ci and C 2 have equal magnitudes when the con- 
straint is obeyed. Again, constraint force information is obtained by inspecting a constraint 
equation in vector form rather than by examining a collection of scalar generalized constraint 
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forces, and the relationship of the multiplier to the constraint forces is completely evident. 
Another advantage to expressing Eqs. (53)-(55) in vector form is that they apply in the 
general three-dimensional case, as do Eqs. (50)-(52), and Eqs. (25)-(27). 

The second example in Ref. [13] involves two particles whose velocities are to remain 
equal in magnitude; however, an additional configuration constraint is imposed on Pi and 
P 2 as they are connected by a rod of fixed length. Zekovich observes the velocities are made 
equal in magnitude by placing a blade at the rod’s midpoint and making the edge parallel 
to the rod. An expression having the same form as Eq. (53), written entirely with scalars, is 
offered as a nonlinear nonholonomic constraint equation. As is the case with Zekovich’s first 
example, the nonlinearity is contrived and it can easily be shown that a linear nonholonomic 
constraint equation describes the constraint dictated by the blade. The diagram on the right 
side of Fig. 2b in Ref. [13] shows a constraint force in the direction of ;V v P| and the other 
constraint force in the direction opposite to ,Y v Pa ; this result can be made to agree with 
Eqs. (55) by renaming the two particles. 

Jankowski has developed an approach for dealing with constraint equations that are not 
necessarily linear in acceleration. A procedure is set forth in Ref. [21] for forming dynami- 
cal equations of motion in which Lagrange multipliers do appear, and then the multipliers 
are eliminated by employing an orthogonal complement matrix to obtain a reduced set of 
equations. The paper concludes with an example involving a single particle P. It is readily 
demonstrated that Eqs. (7) and (8) can be used to obtain the results reported in Ref. [21] 
when the magnitude of the velocity N v p of P in N must have a prescribed time history; 
that is, iv v p • N v p — v(t) 2 = 0. Moreover, inspection of this constraint equation at the 
acceleration level indicates the constraint force applied to P is in the direction of iV v p , and 
Jankowski reaches the same conclusion. However, Eqs. (7) and (8) are not applicable to the 
subsequent case in which the magnitude of the acceleration N a p of P in N is a prescribed 
function of the time t, N a p • N a p — a(t) 2 = 0 

6. Appell’s Particle 

As mentioned earlier, the literature contains ample discussion of an example proposed by 
Appcll in which a single particle must move in a uniform gravitational held so as to satisfy 
an inherently nonlinear nonholonomic constraint equation. A constraint force is identified 
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in connection with this example, and a final brief demonstration of the use of Eqs. (7) and 
(8) shows that they lead to results obtained by Smith 3 and Van Dooren (Ref. [23]). 


6.1. Identification of Constraint Force 

Three motion variables u\, m 2 , and M3 are introduced such that the velocity A v p in a 
Newtonian reference frame N of a particle P is written as 

N v p = Mini + m 2 n 2 + M 3 n 3 (56) 

where fR, n 2 , and n 3 are a right-handed set of mutually perpendicular unit vectors fixed in 
N. Appell’s restriction on the velocity of P is often expressed by the relationship 

u 3 2 = a 2 { Ul 2 + u 2 2 ) (57) 

where a is a constant, ft is pointed out by Smith that the relationship describes a requirement 
for the angle 7 between A 'v p and n 3 , the vertical direction, to remain constant. In fact, the 
constant a is cos 7/ sin 7. The nonlinear nonholonomic constraint equation is differentiated 
with respect to time to bring it to the acceleration level 

2 m 3 m 3 = 2 a 2 (uiiii + u 2 u 2 ) (58) 


where it is linear in U \ , m 2 , and M3; it can be rewritten as 


TV P 


n 3 (Mi 

«3 


TV „ P 


ni + u 2 N sl p • n 2 ) = 


n „ p 


n 3 (Mini + M 2 n 2 ) 

«3 


= 0 (59) 


where 


N n P 


is the acceleration of P in N. Inspection of this equation according to Eqs. 


(5) and (6) indicates that a constraint force C must be applied to P such that the force is 
parallel to the vector within the square brackets; that is, 


C = A 


n 3 (Mini + M 2 n 2 ) 

u 3 


(60) 


This result is in agreement with what is presented by Smith, who shows that C 


N, r P 


v p = 0 


when 


TV,. P 


obeys the constraint. The advantage of inspecting Eq. (59) and immediately 


obtaining the vector form in Eq. (60) is readily apparent; the result is hardly obvious. 


3 C. V. Smith, Jr., “Comments on Geometric Constraints, Virtual Displacements, and Ideal Constraint 
Forces,” private communication, Sept. 25, 2002. 
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6.2. Constraint Force In Evidence 


If one is interested in obtaining a time history of A, equations of motion can be formed 
by carrying out the instructions contained in Eqs. (7). To begin, inspect Eq. (56) to identify 
the vector coefficients of the motion variables. 


N~ P 

Vi = ni, 


N~ P 

v 2 = no, 


N~ P 

v 3 = n 3 


(61) 


Next, identify the force acting on P even when the constraint is not satisfied. The 
gravitational force acting on P is denoted by f = —mgn 3 where m is the mass of P and the 
constant g represents the gravitational force per unit mass. Three dynamical equations of 
motion obtained with Eqs. (7) can be written in terms of vectors as ,Y v r p • (f +C — m N a p ) = 
n r . • (f + C — m N a p ) = 0 (r = 1,2,3), or in terms of scalars 

mill = —Xa 2 ui/u 3 , m« 2 = —A a 2 u 2 /u 3 , rnu :i = A — mg (62) 


in which case they resemble certain expressions found by Smith. When one substitutes u 3 
obtained from the constraint equation (57), the results are identical to Eqs. (3.7) of Ref. [23], 


mill = —A , m« 2 = —A , mii 3 = A — mg (63) 

\/Ui lj rU2 1 \/Ui lj rU2 1 

The fourth relationship needed to determine the unknowns hi, h 2 , h 3 , and A is provided by 
Eq. (58); when it is solved for h 3 and substitution is performed in the third of Eqs. (63), one 
obtains 


. a . . . a 

A = mg + m—(uiui + w 2 w 2 ) = m 9 i 9 0 

u 3 sjup + u 2 


\a(ui 2 + u 2 2 ) 


= mg — A a 2 


(64) 


i/«i 2 + u 2 2 

where the second step is made with the aid of Eq. (57) together with the first and second of 
Eqs. (63). A solution for A is now at hand, and it can be used as a replacement in the first 
and second of Eqs. (63) to yield 


, m 9 • 2 

A = — — g = mg sm 7 
1 + a z 

gaui g sin 7 cos 7 tq 

1 (1 + a 2 )y/ui 2 + u 2 2 a/ui 2 + u 2 2 

gau 2 g sin 7 cos 7 w 2 

2 (1 + a 2 )i/u! 2 + u 2 2 Vui 2 + u 2 2 

A variation of Appcll’s example is considered in Ref. [24], where a 


(65) 

( 66 ) 
(67) 

1 and the grav- 


itational force acting on P is replaced with a force of arbitrary magnitude and direction, 


f = m(F x h 1 + Fyh 2 + F 2 n 3 ). The results of applying the method of Udwadia and Kalaba, 
reported in Eqs. (36) therein, are seen to be in agreement with Eqs. (62) and (65) here after 
setting F x = 0, F y = 0, and F z = —g. Conversely, what has been done here with this 
example so far can be suitably modified so as to reproduce their Eqs. (36) when the lone 
independent constraint equation is dealt with. 


6.3. Constraint Force Not In Evidence 

If one is not interested in a time history of A, one should appeal directly to Eqs. (8) and 
forego use of Eqs. (7) altogether. The dynamical equations of motion (66) and (67), which 
do not contain A, are then reproduced. After embedding the acceleration level constraint 
equation (58) in N a p , 


N P a(uihi + u 2 u 2 ) „ 

a = Mini + u 2 n 2 H . _ = n 3 


Vui 2 + u 2 2 

the required nonholonomic partial accelerations of P in N are readily identified to be 


( 68 ) 


NxP 

a i 


= ni + 


au\ 


r n 3) 


N~ P 


= n 2 + 


au 2 


r n 3 


(69) 


These vectors are clearly not the same as the nonholonomic partial velocities of P in N 
recorded in Eqs. (61). The two equations of interest are then produced by referring to Eqs. 


( 8 ), 


Nx P 


(f + C — m N & p ) = 


Nx P 


(f — m N & 1 ) = 0 (r = l,2). Although some effort is 


required because the equations are coupled in hi and ?i 2 , Eqs. (66) and (67) are recovered. 
No Jacobian or orthogonal complement matrices are involved in obtaining the results in this 
fashion. 


7. A System Containing a Rigid Body 

In practice, the set of particles belonging to a system consists of certain subsets that make 
up rigid bodies. It is important to be able to deal easily with such a system. What follows is a 
presentation of the essential steps needed to extend the discussion in Sec. 2 to encompass rigid 
bodies. The results allow one to deal with an inherently nonlinear nonholonomic constraint 
equation such as N co A • N co B = 0, where n lo a and N ui B are the angular velocities in an 
inertial reference frame N of two unconnected rigid bodies A and B respectively. One is 
then in a position to identify the directions of the constraint torques that must be applied to 
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A and B in order to keep N co A perpendicular to N u B . It also becomes possible to derive, 
directly, explicit analytical equations that govern the constrained motion of the two bodies 
even though the equations are devoid of multipliers. 

When particles P\ , . . . , Pp make up a rigid body B, the acceleration N a Pi in N of a 
generic particle P % of B can be written in terms of the angular acceleration N ot B of B in 
N, the angular velocity n oj b of B in N, and the acceleration N a B in N of B *, the mass 
center of B, 

N a Pi = N a B * + N cx B x r^+ N u B x { N u B x r ?: ) (i = l,...,/3) (70) 


where ly is the position vector from B* to Pi. Now, N cx B can be expressed uniquely as 


N a B = J2 N “r B Ur + 


N a B 


( 71 ) 


r = 1 


where N a B is called the rth nonholonomic partial angular acceleration of B in N. Substi- 
tution from this relationship and from Eqs. (11) into (70) yields 


N~ B* 


y N a Pi u r + N af i = y N a B * ii r + 1 ' l a t 

r = 1 


r — 1 


+ ( J2 W + N ctf ) x iy + N u B x ( N u B x n) (i = l,...,P) (72) 


\r=l 


from which one obtains 


and 


N~Pi 


a t = 


N ~ B* 

a f 


+ N OL B X Ti + 


N u B 


N~ p i 

& r 


N a B * + N cc B x 


x( iV w B xr i ) (i = 1, . . . ,j3) 
(r = 1, . . . , c; i = l,...,P) 


(73) 


(74) 


The latter relationship is the nonholonomic partial acceleration analog to nonholonomic par- 
tial velocity expressions like Eqs. (4.6.5) and (4.11.16) in Ref. [30] used in the case of simple 
nonholonomic systems to obtain contributions of B to F r and F*. Hence, the contribution 
of B to F r is given by 


(Fr) B ± y N a„ * • R, 

2=1 

P 


£(V + x r, 


N x.B 
r 


2=1 
N~B 


= iV af ! " • R + N ct B • T 


N~ B * \ ^ T) i IV ~ 

i — a r • 2_^ Ri + 

(r = 1 , . . . , c) 


2=1 


P 

y r i x R ' t 

2=1 


(75) 
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where the set of all contact forces and distance forces R, acting on the particles of B is 
equivalent to a force R whose line of action passes through B*, together with a couple whose 
torque is T. The constraint forces and torques that must be applied to B in order to satisfy 
nonlinear nonholonomic constraint equations may be included in R and T, or they may be 
omitted; in either case they will not contribute in aggregate to F r . With a similar exercise 
the contribution of B to F* is found to be 


TO* = -E 


N~ Pi 


m ■ 


N ^ Pi 


i= 1 

P 


-£( 


N~B * i Nx.B 

d, 


i= 1 


- N z B * 

— — a r 


T + ot r XTi) ■ rn 
P P 


N a Pi 


E m ^' a 

i = 1 


N „ Pi N ~ B 
- 


E V P 

Ti X 77li a 1 


_ N~B* 

— a r 


R* + N ct r B • T 


i = 1 

(r = 1 


(76) 


where R* and T* are, respectively, the well-known inertia force and inertia torque for B in 
N, formed for use with Kane’s method. 

The procedure for obtaining directly a minimal set of dynamical equations of motion 
for a complex nonholonomic system is seen to bear a very close resemblance to Kane’s 
method for simple nonholonomic systems, the only difference being that one uses n Sl b * and 
N 6l b (r = l,...,c), vectors that are distinct from the familiar vectors Nj v b * and n £j b 
(r = 1 

One may be interested in the constraint forces acting on a rigid body, and therefore form 
equations of motion according to Eqs. (7). In that event it becomes desirable to adapt the 
process of inspecting a constraint equation written at the acceleration level so that one may 
identify the direction of a constraint force and the point to which it is applied, together with 
the direction of a constraint torque and the body upon which it is exerted. 

In a constraint equation having the form of (5), the terms associated with P{, , Pq can 
be rewritten: 

XEa Pi -W is + Z s 
2—1 

P 

= J2[ N a Q + N ct B xr,+ N u B x ( N u B x r*)] • W is + Z s 
2—1 
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( 77 ) 


= ».« • £ W tt + • |>« x W* + £[»«* x 

2 = 1 2 = 1 z=l 

= w a°-W s + ff a B .r s + ^ (s = 1, . . . ,£) 


( 


N u> B 


X Tj)] • Wj s + Z s 


where ly is the position vector from a point Q hxed in B to Pi (i — 1, . . . , (3). The point 
Q need not be the mass center of B. As discussed in connection with Eqs. (5) and (6), 
the appearance of the vector W !S in Eqs. (77) requires the application of a constraint force 
C i S = to Pi . After selecting the line of action of W is such that it passes through P u 

and defining the resultants 

p p 

W s = ^w is , C S = ^C JS (s = (78) 

2=1 2=1 

the set of forces Ci S) . . . ,Cp s applied to B is regarded as equivalent to a single force C s 
whose line of action passes through Q , together with a couple whose torque is T s . The 
resultant C s is given by 

P P 

c . = E C« = E A * W - = (s = 1, . . . , £) (79) 

2=1 2=1 

and the torque T s is equal to the moment of Ci s , . . . , Cp s about Q, 


P P 

T s = ^ p x C is = q x A s W is = A s t s (s = 1,...,£) (80) 

2=1 2=1 

where t s is the moment of Wi s , . . . , W p s about Q, 

p 

r a = EriX W ia ( S = l,...,£) (81) 

2=1 

One can therefore inspect a constraint equation written at the acceleration level and 
conclude that the appearance of the dot product N a Q • W s requires that B is subject to a 
constraint force C s = A S W S applied to Q, and the appearance of the dot product N ol b • t s 
means B must be acted upon by a couple whose constraint torque is T s = A s t s (s — 1, . . . 

The contribution of B to Eqs. (7) is thus represented by 


(K)b = 


N~B+ 


R* + n Cj b 


(Fr)B = 


- N v Q 


R + 


N u B 


(r = 1, 


,p) (82) 


where the set of all contact forces and distance forces R,; acting on the particles of B is 
equivalent to a force R whose line of action passes through Q, together with a couple whose 
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torque is T. All constraint forces C s applied to 0 are included in the resultant R. and all 
constraint torques T s exerted on B are included in T. The vectors N vft and N are (Ref. 
[30]), respectively, the rth nonholonomic partial velocity of Q in N and the rth nonholonomic 
partial angular velocity of B in N. If the system S to which B belongs is not subject to 
motion constraints described by equations that are inherently nonlinear in velocity (£ — 0), 
then S' is a simple nonholonomic system and Eqs. (82) become precisely the relationships 
provided in Ref. [30] for such a system. If all nonlinearities in the nonholonomic constraint 
equations are contrived, then S is in fact a simple nonholonomic system and should be 
treated as such. 

8. Conclusions 

In dealing with motion constraints that are expressed at the velocity level with relation- 
ships that are nonlinear in velocity, there is a distinction to be made between nonholonomic 
constraint equations in which the nonlinearity is inherent, and those in which the nonlinear- 
ity is contrived. Methods are proposed in this paper for dealing with equations of the former 
type. 

Certain forces and torques are required to ensure satisfaction of nonholonomic constraint 
equations that are inherently nonlinear in velocity. One may be interested in expressing 
these constraint forces and torques in vector form so that their directions are known, and 
there may also be interest in knowing the specific points at which the constraint forces must 
be applied or the particular bodies upon which the constraint torques are to be exerted. 
In that case, one may write constraint equations at the acceleration level in vector form, 
in terms of dot products of vectors, and determine the desired information by the simple 
process of inspection. Such information is not available from any of the methods found in the 
existing literature, where constraint equations are invariably expressed in scalar or matrix 
form. The methodology presented herein provides the information readily and stands as one 
of the paper’s main contributions. As demonstrated here by several examples, this method 
is especially advantageous in cases where the required direction of a constraint force is not 
otherwise obvious. 

When one wishes to construct equations containing evidence of constraint forces and 
constraint torques, solution of which yields time histories of those forces and torques, one 
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forms equations of motion with Kane’s method as though a simple nonholonomic system is 
involved. The vector expressions for constraint forces and torques obtained by inspection 
are included with the vector expressions for the usual applied forces and torques. In this 
way, generalized constraint forces are obtained by using a fundamental definition involving 
dot products of vectors, rather than by forming the product of a Jacobian matrix and an 
array of multipliers as recommended in the current literature. The multipliers introduced 
here bear a clear relationship to constraint force and torque vectors, whereas this is not the 
case with other methods. 

On the other hand, when the constraint forces and torques in question are not of in- 
terest, one may form equations of motion that do not involve those forces and torques in 
any way. Such equations can be constructed directly, in explicit analytical form, without 
first formulating equations that do contain evidence of the constraint forces and torques. 
This is accomplished by employing vectors known as nonholonomic partial accelerations and 
nonholonomic partial angular accelerations; these vectors are distinct from the well-known 
nonholonomic partial velocities and nonholonomic partial angular velocities used to form 
Kane’s equations for simple nonholonomic systems, and they are obtained with the same 
simple process of inspection. The use of an orthogonal complement matrix is required when 
one employs existing extensions made to Kane’s method for the purpose of dealing with 
nonlinear nonholonomic constraint equations. Construction of minimal equations of motion 
without resorting to an orthogonal complement represents a significant advantage over such 
approaches, and constitutes another major contribution of the paper. 
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